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Abstract 

The theory of quantum computation can be constructed from the 
abstract study of anyonic systems. In mathematical terms, these are 
unitary topological modular functors. They underlie the Jones poly- 
nomial and arise in Witten-Chern-Simons theory. The braiding and 
fusion of anyonic excitations in quantum Hall electron liquids and 
2D-magnets are modeled by modular functors, opening a new possi- 
bility for the realization of quantum computers. The chief advantage 
of anyonic computation would be physical error correction: An error 
rate scaling like e~ ai , where I is a length scale, and a is some posi- 
tive constant. In contrast, the "presumptive" qubit-model of quantum 
computation, which repairs errors combinatorically, requires a fantas- 
tically low initial error rate (about 10 -4 ) before computation can be 
stabilized. 

Quantum computation is a catch-all for several models of computation 
based on a theoretical ability to manufacture, manipulate and measure quan- 
tum states. In this context, there are three areas where remarkable algo- 
rithms have been found: searching a data base (15), abelian groups (factor- 
ing and discrete logarithm) (19, 27), and simulating physical systems (5, 21). 
To this list we may add a fourth class of algorithms which yield approximate, 



but rapid, evaluations of many quantum invariants of three dimensional man- 
ifolds, e.g., the absolute value of the Jones polynomial of a link L at certain 
roots of unity: ^(e^)]. This seeming curiosity is actually the tip of an 
iceberg which links quantum computation both to low dimensional topology 
and the theory of anyons; the motion of anyons in a two dimensional sys- 
tem defines a braid in 2 + 1 dimension. This iceberg is a model of quantum 
computation based on topological, rather than local, degrees of freedom. 

The class of functions, BQP (functions computable with bounded er- 
ror, given quantum resources, in polynomial time), has been defined in three 
distinct but equivalent ways: via quantum Turing machines (2), quantum 
circuits (3,6), and modular functors (7,8). The last is the subject of this 
article. We may now propose a "thesis" in the spirit of Alonzo Church: all 
"reasonable" computational models which add the resources of quantum me- 
chanics (or quantum field theory) to classical computation yield (efficiently) 
inter-simulable classes: there is one quantum theory of computation. (But 
alas, we are not so sure of our thesis at Planck scale energies. Who is to 
say that all the observables there must even be computable functions in the 
sense of Turing?) 

The case for quantum computation rests on three pillars: inevitability — 
Moore's law suggests we will soon be doing it whether we want to or not, 
desirability — the above mentioned algorithms, and finally feasibility — which 
in the past has been argued from combinatorial fault tolerant protocols. (It 
is a quirk of human optimism that these petitions are usually pleaded inde- 
pendently: e.g., "feasibility" is not seen as a precondition for "inevitability".) 

Focusing on feasibility, we present a model of quantum computation in 
which information is stored and manipulated in "topological degrees of free- 
dom" rather than in localized degrees. The usual candidates for storing in- 
formation are either localized in space (e.g., an electron or a nuclear spin) or 
localized in momentum (e.g., photon polarization.) Almost by definition (see 
"code subspace" below) a topological degree of freedom is protected from lo- 
cal errors. In the presence of perturbation (of the system Hamiltonian) this 
protection will not be perfect, but physical arguments suggest undesirable 
tunneling amplitudes between orthogonal ground states will scale like e~ al , 
where / is a length scale for the system, e.g., the minimum separation main- 
tained between point-like anyonic excitations. We will return to this crucial 
point. 

But let us take a step backward and discuss the standard quantum circuit 
model and the presumptive path toward its physical implementation. To 



2 



specify a quantum circuit T, we begin with a tensor product Cf (g> • • • ® C 2 
of n copies of C 2 , called qubits. Physically, this models a system of n non- 
interacting spin=| particles. The circuit then consists of a sequence of K 
"gates" Uk,l < k < K, applied to individual or paired tensor factors. A gate 
is some coherent interaction; mathematically it is a unitary transformation 
on either Cf or C 2 ® C 2 , 1 < i, j < n, the identity on all remaining factors. 
The gates are taken from a fixed finite library of unitary 2x2 and 4x4 
matrices (with respect to a fixed basis {|0), |1)} for each C 2 factor) and must 
obey the surprisingly mild condition, called "universality", that the set of 
possible gate applications generates the unitary group U(2 n ) densely (up 
to a physically irrelevant overall phase.) Popular choices include a relative 
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on one and two "particles" , respectively. It is known that beyond the density 
requirement the particular choice of gates is not too important. Let Wr = 
\XiL\ Ui denote the operator effected by the circuit T. It is important for 
the fault tolerance theory that many gates can be applied simultaneously (to 
different qubits) without affecting the output of the circuit W r (|0)<8>- • -<8>|0)). 

Formally, information is extracted from the output by measuring the first 
qubit. The probability of observing |1) is given according to the usual axioms 
of quantum mechanics as: 

P (r) = (OlWjlWrlO), (1) 

where Hi is the projection to |1), ^ 1 ) ' a PP^ e< ^ ^° ^ e ^ rs ^ ( l u ^^- 
decision problem, such as finding the k-th binary digit of the largest prime 
factor of an integer x, can be modeled by a function F on binary strings, 
F : {0, 1}* — > {0, 1}; in our example, the input string would encode (x, k). 
We say F belongs to BQP if there is a classical polynomial-time (in string 
length) algorithm for specifying a "quantum circuit" T(y) (in the example 
y = (x, k)) which satisfies: 

p(r(y))>f if F(y) = l and 
p(r(y))<§ if F(y) = 0. 

This definition suggests that one needs to make an individual quantum 
circuit to solve each instance of a computational problem. However, it is 
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possible to construct a single circuit to solve any instance of a given BQP 
problem of bounded size, e.g., factor integers with < 1000 digits. Moreover, 
by (19) there is a universal circuit, univ.(n,k), which simulates all circuits of 
size k on n qubits: 

w univ . (n , k) (|o • • • o) <g> |r>) = w T \o ■ ■ ■ o) ® |r>. 

Yet another definition allows one to do measurements in the middle of 
computation and choose the next unitary gate depending on the measurement 
outcome. This choice will generally involve some classical Boolean gates. 
This scheme is called adaptive quantum computation. In certain cases, it can 
squeeze general BQP computation out of a gate set which is not universal. 

Implementation of a quantum computer 

It is not possible to realize a unitary gate precisely, and even when we 
do not intend to apply a gate, the environment will interact with the qubits 
causing decoherence. Both imprecision and decoherence can be considered 
in a unified way as "errors" which can be quantified by a fidelity distance 
(17) or a super-operator norm (19). A crucial step in the theory of quantum 
computing has been the discovery of error- correcting quantum codes (28) 
and fault-tolerant quantum computation (25,29). These techniques cope 
with sufficiently small errors. However, the error magnitude must be smaller 
than some constant (called an accuracy threshold) for these methods to work. 
According to rather optimistic estimates, this constant lies between 10~ 5 and 
10 -3 , beyond the reach of current technologies. 

The presumptive path toward a quantum computer includes these steps: 
(1) build physical qubits and physical gates; (2) minimize error level down 
below the threshold; (3) implement decoherence-protected logical qubits and 
fault-tolerant logical gates (one logical qubit is realized by several qubits 
using an error- correcting code.) As a counterpoint, the theme of this arti- 
cle is that implementing physical qubits might be redundant. Indeed, one 
can "encode" a logical qubit into a physical system which is not quite a 
set of qubits or even well separated into subsystems. Such a system must 
have macroscopic quantum degrees of freedom which would be decoherence- 
protected. A super-conducting phase or anything related to a local order 
parameter will not work for if a degree of freedom is accessible by local mea- 
surement, it is vulnerable to local error. However, there is another type of 
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macroscopic quantum degree of freedom. It is related to topology and arises 
in collective electronic systems, e.g. the fractional quantum Hall effect (13) 
and most recently in 2D cuprate superconductors above T c (12, 23). 

Though much studied since the mid-1980's the connection between frac- 
tional quantum Hall effect and quantum computation has only recently been 
realized (11 , 20). It was shown by (4) that the ground state of the v — | 
electron liquid on the torus is 3-fold degenerate. This follows from the fact 
that excitations in this system are abelian anyons: moving one excitation 
around another multiplies the state vector by a phase factor e %< ^ (in this case 
= ^y). The process of creating a particle-antiparticle pair, moving one of 
the particles around the torus, and annihilating it specifies a unitary opera- 
tor on the ground state. By moving the particle in two different directions, 
one obtains two different unitary operators A\ and A 2 with the commuta- 
tion relation A 1 A 2 A^ 1 A^ 1 = e 1 ^, implying a ground state degeneracy. This 
argument is very robust and only requires the existence of an energy gap 
or, equivalently, finite correlation length l . Indeed, the degeneracy is lifted 
only by spontaneous tunneling of virtual excitations around the torus. The 

resulting energy splitting scales as e ; o , where I is the size of the system. 
Interaction with the environment does not change this conclusion, although 
thermal noise can create actual excitation pairs rather than virtual ones. 

Both the ground state degeneracy on the torus and the existence of anyons 
are manifestations of somewhat mysterious topological properties of the v = 
| electron liquid itself. Anyons can be regarded as "topological defects" 
similar to Abrikosov vortices but without any local order parameter. The 
presence of a particle enclosed by a loop on the plane can be detected by 
holonomy — moving another particle around the loop. 

At v — |, the electron liquid on the torus could be used as a logical 
"qutrit" (generalized qubit with 3 states). Unfortunately, this will hardly 
work in practice. Besides the obvious problem with implementing the torus 
topology, there is no known way to measure this logical qutrit or prepare it 
in a pure state. 

A more flexible and controllable way of storing quantum information is 
based on nonabelian anyons. These are believed to exist in the v — | frac- 
tional quantum Hall state. According to the theory (22, 24), there should be 
charge | anyonic particles and some other excitations. The quantum state 
of the system with 2n charge \ particles on the plane is 2 n ~ ^degenerate. 
The degeneracy is gradually lifted as two particles come close to each other. 
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More precisely, the 2 n ~ ^dimensional Hilbert space H n splits into two 2 n ~ 2 - 
dimensional subspaces. They correspond to two different types of charge | 
particles which can result from fusion of charge | particles. Thus observing 
the fusion outcome effects a measurement on the Hilbert space H n . This 
model supports adaptive quantum computation when surfaces of high genus 
are included in the theory and admits a combinatorial description (1) appar- 
ently in the same universal class as the fractional quantum Hall fluid. 

Beyond this, a discrete family of quantum Hall models exists (26) based 
on k + 1-fold hard-core interaction between electrons in a fixed Laudau level 
which appears to represent the same universality class as Witten-Chern- 
Simons theory for SU(2) at level=/c (33). Anyons in these models behave 
as topological defects of a geometric construction (10) and their braiding 
matrices have been shown to be universal (8, 9) for k > 3, k ^ 4. 

Code spaces and quantum media 

Even after the particle types and positions of anyons are specified, there 
is an exponentially large (but finite dimensional) Hilbert space describing 
topological degrees of freedom. In combinatorial models, this Hilbert space 
becomes a "code subspace" W of a larger "quantum media" Y . Thus a 
fundamental concept of cryptography is transplanted into physics. Let V be 
a finite dimensional complex vector space, perhaps C 2 , and Y = V <g> • • • <8> V 
an ra-fold tensor product (where n is typically quite large.) Let W C Y be a 
linear subspace. We call W C Y k-code if and only if: 

w n -^f W 

is multiplication by a scalar whenever O is a /c-local operator (an arbitrary 
linear map on any fc-tensor factors of Y and the identity on the remaining 
n — k factors ) and is the orthogonal projection onto W . We think of 
such spaces as resisting local alteration and in the usual interpretation of Y 
as the Hilbert space of n particles, quantum information stored in W will be 
relatively secure. It is a theorem (14) that the quantum information in W 
cannot be degraded by errors operating on fewer than | of the n particles. 

Let us define a (discrete) quantum medium to be a tensor product Y = 
®iVi as above, where now the set of indices {i} consists of points distributed 
on a geometric surface T, together with a local Hamiltonian H = J^Hi 
(each Hi is a Hermitian operator defined only on those tensor factors whose 
index is within e)0 of the i-th point in the geometry of the surface.) Local 
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Hamiltonians H have been found (10, 20) with highly d-degenerate ground 
states corresponding to modular functors (31,32), (and thus braid group 
representation (16) and link polynomials (18).) In these cases, the ground 
state G of H will be fc-code for k ~ injectivity radius of T ~ V area T. The 
topological degrees of freedom referred to above reside in G. But we do not 
attempt a precise mathematical definition of topological degrees of freedom 
since we would like it to extend beyond discrete system, e.g., to fractional 
quantum Hall ground states. 

In the case when T is a disk D with punctures — physically anyonic 
excitations — a sequence of local modifications to H (see (10)) effects a dis- 
crete 1-parameter family H t of Hamiltonians, where the ground states G tl and 
G ti+1 at consecutive time steps differ by a |_§ J -local operator Oi, Oi(G tl ) = 
Gf.^, . Note that if G+. and Gf... are both /c-code that for O; as above, 
Oi\c t . must be unique up to a scalar (for a distinct 0\, consider the re- 
striction of the /c-local operator o 0\ projected to G ti , n Gt . o o 0[). 
This uniqueness property forces this discrete-cryptographic transport of code 
spaces to coincide (up to phase) with the differential geometric notion of adia- 
batic transport — integration of the canonical connection in the "tautological" 
bundle of <i-planes in Y = C 2 ™. If the 1-parameter family is a closed loop, 
a projective representation of the braid group on U(code) is obtained. By 
choosing H, these can be engineered to be precisely the Hecke algebra repre- 
sentations {p\} associated to the Jones polynomial. From H one can build 
a concrete model of quantum computation; the model and its connection to 
the Jones polynomial are described below. Although the H found in (10) is 
enormously cumbersome, it appears to share the universality class "Witten- 
Chern-Simons theory of SU(2) at level 3" with a simple 4-body Hamiltonian 
(26), which has been proposed as a model for the fractional quantum Hall 
plateaus at v — ^ and y. 

An anyonic model for quantum computation 

A family of unitary representations of all mapping class and braid groups 
with certain compatibility properties under fusion is known as a unitary 
topological modular functor. To define our model, we take only the planar 
surface portion of the simplest universal modular functor, Witten-Chern- 
Simons SU(2) modular functor at level 3, and this reduces to the Jones 
representations of braids {p\}. For an appropriate family of local Hamiltoni- 
ans H t , these representations describe the adiabatic transport of the lowest 
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energy states with n anyonic excitation pairs, these states form a subspace 
W of dimension Fibonacci(2n) — as the 2n anyonic excitations are "braided" 
around each other in 2+1 dimensional space time. In this theory, (we de- 
note it CS5 because of its link to the Jones polynomial at a fifth root of 
unity), there are 4 label types 0,1,2,3 corresponding to the complete list of 
irreducible representations of the quantum group SU(2) 5 of dimensions 1,2,3 
and 4 (or equivalently the irreducible positive energy loop group representa- 
tions of £SU(2) at level 3). We initialize our system on the disk D in a known 
state by pulling anyonic pairs out of the vacuum. This theory is self-dual 
so the two partners have identical types. Pairs are kept or returned to the 
vacuum according to the results of local holonomic measurement. Finally we 
have a known initial state in the disk with 2n punctures where each puncture 
has label=l and dD has label=0. We assume n is even and group the punc- 
tures into n/2 batches of 4 punctures each. Similar to (8), each batch B is 
used to encode one qubit = C 2 : the basis {|0), |1)} is mapped into the type 
(0 or 2) of the 2-fold composite particle (round circle), which would result 
from fusing a (fixed) pair of the type 1 particles within B. Initially, both the 
double and 4-fold composites (ovals) have type 0. By maintaining, at least 
approximately, this condition on the ovals after the braiding is complete, we 
define a "computational summand" of the modular functor: it is spanned by 
n-bit strings of 0's and 2's residing on the 2-fold composites (round circles). 




Figure 1 



Now as in the quantum circuit model, a classical poly-time algorithm 
looks at the problem instance (F, y) and builds a sequence of "gates" , but 
this time the gates are braid generators (right half twist between adjacent 
anyons) , 1 < i < 2n — 1, and a powerful approximation theorem (19 , 30) 
is used to select the braid sequence which approximates the more traditional 
quantum circuit solving (F,y). So the topological model may be described 
as: 

1. Initialization of a known state in the modular functor. 

2. Classical computation of braid b effecting a desired unitary transfor- 
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mation X of the computational subspace of the modular fucntor. 

3. Adiabatic implementation of the braid by (somehow) moving the anyons 
in D to draw b in space-time. (Here we keep the anyons separated by 
a scale I.) 

4. Application of a projection operator II to measure the type (0 or 2) of 
the "left most" composite particle (as seen in Figure 1.). 

The last step is the direct analogue of measuring the first qubit in the 
quantum circuit model and the same formula (1) applies: the probability of 
observing type (the null particle) is: 



where the 0's on the right hand side represent our carefully prepared initial 
state with 2n type 1 excitations. 

To close the topological discussion, we note that the previous formula 
can be translated (using the S matrix) into a plat closure (See Figure 2) of 
a braid L = plat(6~ 1 7&), where 7 is a small loop inserted to measure the 
left-most qubit, and now the outcome of the quantum circuit calculation, 
prob(0), becomes a Jones evaluation 



prob(0) = (OIX+ITXIO), 




minima (L) 



Figure 2 
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where [2] 5 = 1 5 ; c, m, and iu are the number of components, number of 
local minima, and writhe of L respectively; and a = e 7 ™/ 10 . Details of this 
calculation will be posted at url: www.tqc.iu.edu. 

The braiding disturbs and reforms composite particle types with sufficient 
subtlety to effect universal computation. To reduce this model to engineer- 
ing, very significant obstacles must be overcome: stable quantum media must 
be maintained in a suitable phase, e.g., CS5; excitations must be readily ma- 
nipulated, and electrical neutral particle types and 2 must be distinguished, 
presumably by holonomy experiments. Although these challenges are daunt- 
ing, they are, perhaps, less difficult than a head-on assault on the accuracy 
threshold in the quantum circuit model. 
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